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Abstract 

The aim of this work is to describe the electronic properties of graphene in a constant magnetic 
field in the long wavelength approximation with random binary disorder, by solving the Soven 
equation self-consistently. Density of state contributions for different valleys in each sublattice sites 
are obtained for different values of magnetic field strength showing remarkable differences between 
K and K' valleys. A band gap is obtained by an asymmetric on-site impurity concentration and 
the graphene electrons acquire an anomalous magnetic moment, which is opposite in different 
valleys, which depend highly in the interplay between the impurity band, the band edges and the 
broadening of the Landau levels. In turn, magnetization as a function of B for different on-site 
random impurities is computed showing that by decreasing the on-site impurity energy values, 
maximum magnetization is shifted towards higher values of B which can be used to create and 
manipulate polarized valley currents. Finally, conductivity and local vertex function are obtained 
as a function of energy showing that scattering contributions from A and B sublattices differ 
significantly. Effective medium local two-irreducible vertex is computed showing that scattering 
from sublattice A and B do not contribute equally, which can be related to weak anti-localization. 
From these results, it could be possible to explore how the valley pseudospin can be used to create 
polarized currents by populating asymmetrically the sublattice sites, where the population can 
be tuned with the applied magnetic field strength. 


1 Introduction 

Graphene is a two-dimensional allotrope of carbon which has become one of the most significant topics 
in solid state physics due to the large number of applications mmm , m ed- The carbon atoms 
form a honey-conrb lattice made of two interpenetrating triangular sublattices, A and B. A special 
feature of the graphene band structure is the linear dispersion at the Dirac points which are dictated 
by the n and n' bands that form conical valleys touching at the two independent high symmetry 
points at the corner of the Brillouin zone, the so called valley pseudospin [6]. In the absence of 
defects, electrons near these symmetry points behave as massless relativistic Dirac fermions with an 
effective Dirac-Weyl Hamiltonian [Jj. When a magnetic field is applied perpendicular to the graphene 
sheet, a discretization of the energy levels is obtained, the so called Landau levels [Jj. These quantized 
energy levels still appear also for relativistic electrons, just their dependence on field and quantization 
parameter is different. In a conventional non-relativistic electron gas, Landau quantization produces 
equidistant energy levels, which is due to the parabolic dispersion law of free electrons. In graphene, 
the electrons have relativistic dispersion law, which strongly modifies the Landau quantization of the 
energy and the position of the levels. In particular, these levels are not equidistant as occurs in a 
conventional non-relativistic electron gas in a magnetic field. This large gap allows one to observe the 
quantum Hall effect in graphene, even at room temperature [8]. In turn, the valley pseudospin can 

* email: jsardenghi@gmail.com, fax number: +54-291-4595142 


1 



be used to create polarized currents in a similar way as the real electron spin is used in spintronics. 
By applying a local gate voltage to a quantum nano-contact in a graphene nanoribbon, valley filter 
can be obtained. Due to the finite size of the graphene sheet, the transvesal momentum component 
are quantized. Those states with a definite group velocity, depending on the polarized gate voltage 
sign applied, can go through the quantum nano-contact with the result of a largest population in one 
of the valleys. Altough the theoretical approach to distinguish carriers in the two valleys are known, 
the experimental procedure has become an attracting literature (see [9], [TO], [XT] and [12)- This 
valley pseudopin carries real magnetic moment, like real spin, supported by chiral orbital current, 
which can rise Zeeman splitting and Pauli paramagnetism in graphene m This magnetic moment 
can be used to couple it with a perpendicular magnetic field, which can be about 30 times the Bohr 
magneton for low energy electrons [9]. However, taking into account the total contributions of the 
valley pseudospin, there is no net polarization by the magnetic field because the contribution of the 
Zeeman term to the Hamiltonian has opposite sign at the K and the K' points, which implies that 
the direction of the valley pseudospin polarization induced by the magnetic held cancel each other 
M- In this case, the valley degeneracy is protected by spatial inversion symmetry m- Besides 
the usual spatial translation and rotation symmetries, pristine graphene Hamiltonian is invariant 
under a large number of symmetries in the isospin spaces [16] . By introducing the most general 

disorder Hamiltonian Hdi S = pseudospin symmetries can be broken depending on which 

ij 

elements of Vy are non-zero ([T7j and [15])- In particular, local diagonal disorder can be obtained 
experimentally by irradiation m, where carbon atoms are extracted from the graphene plane and 
adatoms or adsorbed species attach to the graphene plane [20J . or in which some carbon atoms are 
chemically substituted for other elements. These adsorbed particles can induce a local potential at 
the sites where they couple to the carbon lattice and change the on-site energy in the Hamiltonian. 
For weak random substitutional disorder potential, intravalley mixing within either the K and K' 
valleys is possible, but not intervalley mixing. In the case of a random binary alloy disorder, a gap 
can be open in the Fermi level due to the broken C z chiral symmetry m- The peculiar behavior of 
the n = 0 Landau level in graphene, where its amplitude is nonzero only in one of the sublattices, 
namely, at B sites for K valley and A sites for K' valley, combined with the random diagonal disorder 
can enhance the asymmetry in the valley amplitudes, which is related to the geometric nature of the 
Bloch band and its relation with Berry phase (see [21], [22] . [ 231 . [23] and [21] )• This asymmetry can 
be detected by a population difference in the two valleys as a signal of orbital magnetization (see [Til ) 
and in turn, this asymmetry can be enhanced by an asymmetry in the relative density of states for 
both valley pseudospins. In this sense, the aim of this work is to study the behavior of the valley 
pseudospin under a constant magnetic field and random binary disorder. In particular, self-energies 
and density of states for both sublattices and both valleys will be obtained by applying coherent 
potential approximation (CPA) (see |2B], [27] and [28]). In this method, the system is replaced by an 
effective medium with a complex self-energy that replace the disorder random potential. The value 
of the effective medium self-energy can be obtained by demanding that the average scattering of an 
electron by the surrounding medium is zero. Improvement of the method for better understanding of 
non local disorder can be done by replacing the single impurity site by a cluster of atoms (see [29j . 
1501 and E5). Due to the unbroken translational symmetry, the embedded cluster method has been 
applied ([32], ESI and E3)- To obey all the imposed criteria for the CPA (see [3S]), the Non Local 
Coherent Potential approximation has been developed m which is based in the approximation of 
continuous lattice functions f(k) by cluster function f(K n ), where K n are points in momentum space 
that satisfy Born-von Karrnan boundary conditions. Without taking into account the improvements 
of the CPA method, we will consider the results for low energy Bloch electrons with self-energies in the 
full Born approximation with nested diagrams included]]] This work will be organized as follow: In 
section II, the magnetic Green function with diagonal on-site energies will be computed for graphene. 
Although this procedure has been studied without diagonal on-site energies, for a self-contained work 
this results will be generalized. In section III, single-site approximation will be applied and a system 
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of coupled Soven equation will be found and solved. The discussion of the results is shown in section 
IV and the principal findings of this paper arc highlighted in the conclusion. 


2 Green function 


For a self-contained lecture of this paper, a brief introduction of the quantum mechanics of graphene 
in a constant magnetic field in the long wavelength approximation will be introduced (see [38]). The 
Hamiltonian in the two inequivalent corners of the Brillouin zones can be put in a compact notation 
as 


Hq X) = V F (Xa x P X + CTyPy ) + 


VA) 

-‘A 

0 


0 

VA) 


(1) 


where A = 1 is for the K valley and A = —1 for the K' valley. A diagonal energy matrix has been 
introduced for further application to single-site approximation (CPA). The quasiparticle momentum 
is p—eA, where e is the electron charge and A is the vectorpotential which in the Landau gauge 
reads A =(— By, 0,0) and v F = 10 e m/s is the Fermi velocity]^ The Hamiltonian of last equation is 
invariant under traslations in the x direction which means that the wave functions can be written 
ikx a coordinate transformation y = hv F k x — eBy can be applied followed 

y , finally, introducing the annihilation and creation operators 


as i />( A ) = e 


by a scale transformation y = 


/ ehB 


-j= [y + J= j and ed = (y — J= j , the Hamiltonian of eq.® reads 


#(A=i) = 


Ai) 


7 a 1 



#(*=-!) = 


A-i) 

J A 
—7 a 


-0 


( 2 ) 


where 7 = v F y/2ehB. The eigenfunctions and eigenvectors for the Hamiltonian of last equation reads 
(see [39] and (IQ]) 


il> 


(A) 

(n,s,fe) 


(r) = e 


ikx C' n , 


V2L^ { 


n,s,k ) 


(0 


where reads 


P(n, a ,k)(S) = 


_( an’ S Vn-l,fc(0( 1 “ <*n,o) 


4 > n,k{£) 


V’jni!fc)(0 = 


being <p n ,k{£,) the wave function of the harmonic oscillatojfl 


0n,fe(C) 

c4r 1,s) 0n-l,*(£)(! - <5n,o) 


(3) 


(4) 




7 T - 1 / 4 

y/2 n n\ 




(5) 


= l sk,L = VA where A is the area of the graphene sheet and 


V (A>«) 


A(s7 - Z ( B X) ) - S] /(S ( A X) - E ^) 2 + *y*r 
27 y/n 


( 6 ) 


2 Spin Zeeman energy are completely neglected because the spin splitting is much smaller than Landau-level separa¬ 
tions. 

3 The factor (1 — <5 n? o) is introduced to discriminate the wave function with n = 0. In this case, only one sublattice 


contributes in both valleys K and K'. 


3 










where s = ±1 is the conduction (valence) band index. The coefficient C n is C n = 


•y/2—<5 n ,o 


l b = sJh/eB is the magnetic length. The eigenvalues of the Hamiltonian reads 

E(s) = + 4 A) - V(4 A) - *4 A) ) 2 + ^ 


and 


(7) 


The degeneracy of each level is g = L 2 B/cf>o where <po = h/e is the quantum of flux. The low energy 
description is only valid as long as the characteristic energy of the excitations is not larger than an 
energy cutoff E n> i < Ec, where Ec = Tivpk a and is a momentum cutoff. We can choose k& in 
such a way to conserve the total number of states in the Brillouin zone, that is, irk\ = (27 t) 2 /Ac, where 
Ac = 3V3a 2 /2 is thea area of the hexagonal lattice (see [41]). Then, using eq.0, E < Ec implies that 

s _— 1 _|_ w here S = yJ 87t/ 3-\/3, then for weak magnetic fields, the cutoff 

tends to infinity and for high magnetic fields, the cutoff tends to zero. Eq. 0 indicates that asymmetry 

in the substitutional impurity energies opens an energy gap = yj (E^ — E^ •*) 2 + 47 2 n. 

Using the spectral representation, the Green function of this system reads 


n < 


2 ehB 


,(A) 


+«A 


( r y,E) = — y / E 


f+OO 


(a^ A,s) ) IA+i—i-lI e »fc(x-x') ±* 
v ' n-i-\ —r—. 


2L n—0 «=±l- , -°° 

A,*(0(1 - 


( 8 ) 


(2 - S n , 0 ) (E - E ( n s) ) 


where for the moment i, j = 0,1 where 0(1) represent the sublattice A(B). We can perform the 
integration in k by using eq. and completing squares 

ii '\ 1 ( v ink ) 2 1 ( y' ink ) 2 — l%(k— y + v UU —^-l ) 2 —k-[(y—y , ) 2 + (x—x 1 ) 2 ]^ — h-Py+y')(x—x / ) 

e lk(x-x ) e -2(TJ- L Bk) e ~2^T_ e 2lJ ) e 4 U ^ 1 A ’ /g\ 

Introducing the following coordinate transformation 

y + ff + i(x — x 1 ) 


q = -Ink 


21 T 


and using the following relation (see eq.(7.377 8 ) of page 804 of [43]) the Green function reads 


£> {r y,E) = 


2 L 


( 10 ) 


( 11 ) 


( 1 -^, 0 ) lx+1 - i ~ jl T( x, i , j) J 

h ( 2 -<w>) " V 


#-<(»-j + V)! r |j-<| 


(n-i+1^4)! n -J+-r- 




where the sum in s has been performed and where 


/(r,r') = 


y-y - i[x- x’) 
21 b 


rj(r,r') is a gauge term that reads 


and 


v(r,r') = ^2 -(y + y)( x ~x') 


P = 


| r — r 


( 12 ) 

(13) 

(14) 
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The Tn X ’ 1 ’^ matrix elements reacQ 


r (|A|,0,0) = r (-|A|,l,l) = ~ 2E ^ + (W + E Q ( 1 X) - Q2 X) )Q ( 1 X) ~ Q^E + Qj X) 

n'y 2 [-E 2 +n'y ‘ 2 + 'E ( x\E-Y, {X) ) + ET, {X) ] 


t( a,i,o) = t( a,o,i) = - W - A(sj A) - Z^>)[E + Q\ a> - S£ ; - S^] 
7 ^[™ 7 2 - E 2 + ?, {X \e - Y, {X) ) + EY, {X) ] 


-(AT 


-(A) 


-(A) 


(A) - 


ji( 1 , 1 , 1 ) _ j,(-l, 0 , 0 ) _ 


-2E+E 


(A) 


-(A) 


n 7 2 - E 2 + T, {X) (E - E^ A) ) + EE {X) 


(15) 


where 

Q[ X) = \[0- + A)(S^) i 5 - (A - l)(S^)i] (16 ) 

in the case that E^ = E^ = 0, eg. (fill) is identical to eq.(ll) and eq.(12) for A = 1 of [i5] . It should 
be pointed out that, as it was shown in eq.(|T]) the application of a magnetic field to graphene allows to 
obtain different wavefunctions for sublattice A and B in each valley. This is different of what occurs 
when B is zero in the effective low-energy description, because in this case the wavefunctions in each 
sublattice site differ only by a phase which depends on the polar angle of the wave vector k. For this, 
the probability amplitudes in both sublattice sites are identical. This crucial difference is of major 
importance when diagonal disorder is introduced because it will increase the availability of states near 
the Fermi energy without disabling the sublattice asymmetry introduced by the magnetic field. 


3 Single-site approximation 


To apply CPA we can introduce impurity potentials in the A and B sublattices in the following form 


H = V F (\(T X P X + (TyPy) + 


Va o \ 

0 V B J 


(17) 


where the potential V A ! B reads 


n a/b 

yA/B = yA/B 


In) (n| 


( 18 ) 


A / d a / R 

and 7 are the on-site impurity energies with a probability distribution P{V i ' ) which is the 
same for both sublattices and pseudospin valleys^ To apply CPA, we can introduce an effective 
Hamiltonian 

H e f = V F (Xa x p x + VyPy) + E (A) (,z) (19) 

where S(z) is a energy-dependent self-energy matrix which reads 


Z w ( z ) 


^\z) 0 \ 

0 S^z) ) 


( 20 ) 


Then, an impurity can be introduced in the effective medium in an specific site ro by defining the 
following Hamiltonian 

H h = H ef + W\r 0 )(r 0 \ (21) 


where W reads 


W = 


V A - E A (z) 
0 


V b -Eb(z) J 


( 22 ) 


4 The coefficients have been introduced to write ea. dllD in a more compact way, although its units are E —1 

which are the units of the Green function. 

5 The random distribution is identical for A and B sublattices. 
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Applying the self-consistency condition for the Green function G e f = ( Gh ), where (...) means config¬ 
urational averaging, we obtain the Soven equation in matrix forn^j 

(w[I-g ef (r o ,r o )W]- 1 ) = 0 (23) 

where g e f is the Green function matrix with magnetic field associated to the effective Hamiltonian of 
eq. dUD, which is identical to eq. m and implies that the effective Green function is that of eq. m, 
which must be evaluated at rg in both coordinates arguments 


9 


W (r 0 ,r 0 ,E) = -fciV (1 ~ 


2 L ^ 

n =0 


(2 — S n fi) 


(24) 


The non-diagonal elements are zero due to the factor f(r,r) = 0 for identical points. This result 
implies that eo. (l23l) contains a system of two coupled self-energies that must be solved self-consistently. 
Applying the configurational averaging by using the probability distribution for random binary alloy 
P(Vi) = cS(V — ei) + (1 — c)8(V — £ 2 ), the system of eq.([23D reads 


c(ei - E^ A) ) 

1 ^ 9aa 

c(e 1 - 

O-S 

-4 A) ) 

1 (A) 

1 9bb 

O - Z ( B } ' 


(l-c)( e 2 -S^) 


1 - 9aa 

f v(^) 

e 2 - iu A 

(1 - c)(e 2 - E^) 

1 (A) 

1 9bb 

f v(^) 

e 2 - 


0 

0 


(25) 


where we have restored the original notation = 9oa an d ffae = 9ii\ where these functions 
depend on E^ and E^ ^ and c is the concentration of e\ impurities that can be located in A or B 
sublattices and e 2 is the on-site impurity energy with probability 1 — c that can be located in the A 
or B sublattice. 


4 Results and discussion 


To apply CPA to ea. CSl) we have to take into account the diagonal elements of the Green function, 
which can be solved exactly 


1 B^ 

Ilcw 


9aa (A A E) = g ( B g ] (r 0 , r 0 , E) = 


4 L 


ip {0 \n A + 1 + C (A) ) - V> (0) (1 + C (A) ) 


^(°)(1 + C W) + ^(0 )( nA + i) _ ^,( 0 )(C-(A) 

+ n A + 1) + 7o 


(26) 


where 70 is the Euler-Mascheroni constant, is the Digamma function and the coefficients A, B 
and C reads 


A {X) (E) = 2q -2 



B (X \E) = 7 -4 
C W {E)= 7 


\EQ[ X) -Q {X) )Q[ X) ~Q {X) E + Q 
2 \-E 2 + S^ A) (E - E^ A) ) + EY, {X) 


(A) 

3 


6 For more details about CPA see EE chapter VIII. 


( 27 ) 
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Figure 1: Total density of states for two different values of applied magnetic field. Blue line for K 
valley and red line for K' valley and ei = —0.5 = —€ 2 - 


In the same way 


where 


9bb{^ r > E ) = 9aa\^ A E ) = 


-2 E + E^ 


AA) 

-‘B 


2 L —E 2 + e[£\e - E ^°) + EE^ 

—-D (a) U (0) (F (a) +n A + 1) -^ (0) (1 + F (a) ) 

4:L 


D W (E) = 7 


-2 


-2E + E 


(A) 


AA) 


F (a) (E) = 7“ 2 -E 2 + E ^ A) (E - E^ } ) + EE ( g } 


A 

(A)> 


(28) 


(29) 


From eg. ESI) and ea. (E51) . the contribution of the conduction (valence) band in the sublattice A ( B ) 
equal the valence (conduction) band in the sublattice B (A), which reflects the fact that, although 
on-site impurities were introduced, the configurational averaging restore the full symmetry of the 
Hamiltonian. A convenient way to treat the coupled system of Soven equations of eq. (ESI) is to 
consider an iterative formula equivalent to the CPA condition in the following way 


Ei A) B (z+l) = E^ B (*) + 


(A) 




^$ B (<)-V + 


C/bbK/bW'CW)) 


(30) 

In order to obtain the value of the self-energies, the iteration described above for both sublattices 


is repeated until 


4 V + 1 )- s a^) 


< s, where £ determines the requested precision of the 


calculation. The total and partial density of states can be computed as 

p{E) = --lmg(E -E{E)) 

7 r 


( 31 ) 
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Figure 2: Projected DOS on A sublattice for two 
different values of magnetic field. Blue line for K 
valley and red line for K' valley and £a = —0.5 = 
—£b 


Figure 3: Projected DOS on B sublattice for two 
different values of magnetic field. Blue line for K 
valley and red line for I\' valley and £a = —0.5 = 
—£b 


In figure [I] and figures [2] and [3] the total density of states and projected DOS in A and B sublattice 
is shown for both pseudospin valleys using a = —0.5 = —£ 2 , c = 0.7 and for B = 10T and B = 40T. 
For the specific concentration chosen, we are considering acceptor impurities. The total DOS shows 
effectively that no broken pseudospin symmetry has been produced by averaging over disorder, but 
figure 03 and [3] shows an asymmetry in the K and K' valleys for the projected density of states below 
the Fermi level in both sublattices. There is an enhancement of the population for A sublattice in the 
K valley with respect to the K' valley and the opposite for B sublattice. This behavior is expected due 
to the asymmetry of the lowest Landau level in both sublattices, where the shifted peak correspond 
to the impurity contribution. The behavior of the density of states near the energy threshold where 
the approximation holds is expected and reflects the linear regime with the peaks corresponding to 
the broadening of Landau levels (see figure 10 of SB)- In turn, the impurities open an energy gap at 
the Dirac point (see [44]) which lies below the Fermi level for both magnetic field strength, although 
for higher B , the impurity band is reduced and shifted away the Fermi level. In the projected DOS 
figures, the energy gap introduced by the diagonal impurities is different for K and K' valley, which is 
a consequence of the self-energy asymmetry results obtained from CPA. In the case that c = 0.5, the 
asymmetry in the gap for K and K' valleys should be of major importance because in this case, the 
impurity band will be located in the Fermi level and for half-filling, a net magnetization will appears^ 
In figure [4] and [5] the net magnetization in both sublattices is shown for different impurity on-site 
energies where 

Ha /B /»o = [Pa/b( e ) - PA /B (E)]dE (32) 

J—Ec 

where /io is the Bohr magneton. As it is expected, both magnetization are equal an opposite, which 
implies that there is no net magnetization when both sublattices are taking into account. In both 
sublattices, the magnetization present a maximum value that decreases with an increasing on-site 
impurity potentials, which implies that it is always possible to maximize magnetization with an 
specific choice of the parameters. The shift of the maximum with the displacement of impurity on-site 
energy is the combination of two physical process: the enhanced asymmetry between K and K' density 
of states near the Fermi level when the strength of the impurities is raised and the spreadening of the 
impurity band due to the weak magnetic field, although when electron many body effects are taking 

7 The c = 0.5 impurity concentration has not been studied in this work because convergence for self-energy in CPA 
is difficult to achieve in this case. 
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Figure 4: Magnetization of A sublattice due to relative Figure 5: Magnetization of B sublattice due to relative 
difference between density of state contributions for difference between density of state contributions for 
the two valleys pseudospin. the two valleys pseudospin. 



Figure 6: Conductivity as a function of energy for different on-site impurities. Blue line £i = e 2 = —0.6, 
green line ei = £2 = —0.8 and red line for = e 2 = — 1. Dashed line for sublattice A scattering 
contribution to conductivity. The magnetic field strength used is B = 10 T. 
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Figure 7: Vertex function as a function of energy for different on-site impurities. Blue line e± = £2 = 
—0.6, green line ei = £2 = —0.8 and red line for = £2 = —1. Thick line for sublattice A scattering 
and dashed line for sublattice B scattering contribution to conductivity. The magnetic field strength 
used is B = 10 T. 


into account, this behavior should be reduced because high density of states available from the K or 
K' band increase the energy caused by spin-alignment which in turn results in a small contribution 
to the exchange energy. 

In figure El the conductivity is plotted against energy using eq.(14) of @5], where Kubo formalism 
is used, for different values of impurity on-site energies. The position of the conductivity peaks 
correspond to the minimum and maximum of the energy band, where the density of states is the largest. 
In dashed line, the contribution to conductivity for sublattice A scattering is taking into account. In 
particular, for higher values of impurity on-site energies, the contribution to total conductivity peaks 
are given by A sublattice scattering (red line), but is not the case for the blue line and green line, 
where there is a contribution from B sublattice scattering for the second peak respectively. In figure 
[3 the local two-particle irreducible vertex A consistent with the coherent potential approximation for 
the self-energy E is plotted against energy, for retarded and advanced Green functions (see ESI, S3. 
H8] and [J5]) with the same parameters as figure [6] The peak of the vertex for A sublattice scattering 
(blue thick line in figure 0) is related to peak near the Fermi level in figure El while the second peak is 
given by the contribution of vertex for B sublattice scattering (blue dashed line in figure [TJ) with its 
respective decrease in the conductivity contribution (see blue dashed line in figure EJ). In turn, from 
figure 0 the vertex function for A sublattice scattering shows a singular behavior where the retarded 
and advanced Green function in the real part of the energy are identical. The same behavior is not 
present for the B sublattice scattering which is related to the different contributions of Landau levels 
for A and B sublattice. The singular behavior of the vertex function implies a critical contribution 
to conductivity that change sign and can be related to weak anti-localization (see [50], [51] , [52], 53] 
and E5). 

These results reflect the fact that, although a simple model has been considered, where the im¬ 
purities are diagonal in the sublattice basis and short-ranged and averaging has been applied, the 
application of a constant magnetic field introduce subtleties in the valley properties of electrons. In 
particular, the different Landau levels contributions to the eigenfunctions for each sublattice increase 
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the asymmetry in the orbital magnetic moment present in electrons in the two Dirac points and for 
particular values of impurity concentrations and magnetic field strength, the net magnetization could 
be larger to create and manipulate polarized valley currents. From the experimental viewpoint, doping 
in graphene can be obtained through electric doping by changing gate voltage (see [S3] and }5B]) or by 
chemical doping, which is discussed as surface transfer doping and subtitutional doping (see m and 
m)- For a general review of the experimental procedure to obtain doping asymmetry in graphene 
see [59] , (60] and the references therein and for the direct experimental determination of the chemical 
bonding of individual impurity atoms see m- In the case that p-type and n-type doping can be 
achieved in a controled way and where the experimental determination of the on-site energies results 
in slightly different values, then it would be possible to tune up population difference in both valleys 
in each sublattice site through the magnetic field strength and valley currents could be obtained by 
applying a gate voltage. 

5 Conclusion 

We have investigated the behavior of the density of states in both sublattices and for the different valley 
pseudospins in graphene with disorder and magnetic field in the low energy effective-mass theory. By 
applying coherent potential approximation for a binary random alloy and considering the contribution 
for both sublattices on the effective self-energies, a coupled system of Soven equations can be obtained 
to be solved numerically. We have shown that a band gap is opened by an asymmetric on-site impurity 
concentration and the graphene electrons acquire an anomalous magnetic moment, which is opposite 
in different valleys, similar to a real spin. The valley contributions to the projected density of states 
for each sublattice are not simetrical and the differences depend highly in the interplay between the 
impurity band, the band edges and the broadening of the Landau levels. In the case that spin- 
orbit interaction and Zeeman effect are considered, the asymmetry between valleys can be enhanced 
by splitting of Landau levels. In turn, conductivity and CPA local two-irreducible vertex has been 
computed showing that scattering from sublattice A and B do not contribute equally, which implies 
that graphene Bloch electrons correlations in the long wavelength approximation in averaged disorder 
and magnetic field develop unusual behavior for critical energy values. 
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